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Abstract. Class sources are objects representing the earliest phase 
of the protostellar evolution. Since they are highly obscured by an ex- 
tended dusty envelope, these objects emit mainly in the far-infrared to 
millimetre wavelength range. The analysis of their spectral energy dis- 
tributions with wide wavelength coverage allows to determine the bolo- 
metric temperature and luminosity. However, a more detailed physical 
interpretation of the internal physical structure of these objects requires 
radiative transfer modelling. We present modelling results of spectral 
energy distributions of a sample of nine Class sources in the Perseus 
and Orion molecular clouds. The SEDs have been simulated using a ra- 
diative transfer code based on the Monte Carlo method. We find that a 
spherically symmetric model for the youngest Class sources allows to 
reproduce the observed SEDs reasonably well. From our modelling we 
derive physical parameters of our sources, such as their mass, density 
distribution, size, etc. We find a density structure of p ~ r~ 2 for the 
collapsing cores at young ages, evolving to p ~ r~ 3 / 2 at later times. 
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1. Introduction 

A Class source is an object representing the earliest phase of 
the protostellar evolution (e.g. Andre et al. 1993). It consists of 
a central protostellar object surrounded by an infalling dusty enve- 
lope and a flattened accretion disk. Many (presumably all) Class 
sources are associated with a bipolar molecular outflow. Because of 
their difficult detection (their photospheres are highly obscured and 
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the objects spend only a short time (~ 10 4 yr) in this evolutionary 
phase), and the lack of knowledge in model protostellar parameters 
(distribution, composition and size of dust, density and temperature 
distribution, etc.), constraining physical properties to Class objects 
is quite difficult. Continuum sub-mm observations of Class sources, 
however, let us detect dust emission of the massive circumstellar 
envelope of these sources. Information from the dust emission, a 
physical model of the extended envelope as a pre-requisite, and the 
implementation of techniques like the blackbody and the envelope 
fitting procedures, open a way to interpret the structure of Class 
sources. 

2. An envelope model 

How can we describe the protostellar emission from the enve- 
lope? We adopt the standard envelope model (Adams 1991) to keep 
the problem as simple as possible. The circular symmetry of the 
observed emission (Rengel et al. 2004a), and the lack of significant 
internal structure justifies the simplicity of a spherical model case. 
If the emission is optically thin, the observed intensity for a sphe- 
rical symmetric protostellar envelope at an impact parameter b is 
given by Eq. 1 assuming single power-laws of the opacity, the radial 
temperature and density distribution ((3, q and p, respectively). 



ro is the outer radius, p is the density, the dust temperature, 
k v the opacity of the dust grains and B v [Td(r)] the Planck function at 
dust temperature TV If the emission is in the Rayleigh- Jeans limit 
and if r » b, Eq. 1 can be approximated to I u (b)/I u (0) = (b/bo)~ m 
(where m is the power-law index of the observed intensity profile, 
and bo is the normalization factor to the peak emission). 

Results from the standard theory provide the first direct insights 
into observable estimations. Nevertheless, the temperature profile 
will diverge from a single power-law q as the envelope becomes op- 
tically thick at the primary wavelengths of energy transport (inner 
portion of the envelope) (e.g. Shirley, Evans & Rawlings 2002). Here 
it becomes necessary to calculate the temperature distribution self- 
consistently by implementation of a radiative transfer code. To cal- 
culate the observables (temperature and density distributions, and 
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Figure 1: Temperature profile (left) and best simulated and observed 
SEDs (right, short dashed and dotted lines, respectively). Details 
regarding the summarized target list at several wavelengths from 
the literature are given in Rengel et al. (2004b). 



intensity map) we implemented a ID spherical symmetric radiative 
transfer code (MC3D, Wolf et al. 1999). The model consists of a 
central heating source, with stellar luminosity and an effective 
temperature T*, embedded within a spherical symmetric dust enve- 
lope of mass M env . The envelope is specified by an outer radius R ou t, 
a density power-law index p, and inner radius r in , which is set by the 
dust destruction radius r su b, adopted to a temperature of 1500 K. 



3. Modeling continuum submillimetre emission 

First constrains to physical properties (spectral indices, masses, 
radial profiles and sizes) of 15 objects are derived from continuum 
SCUBA imaging of six star formation regions in the Perseus and 
Orion molecular cloud complexes (Rengel et al. 2004a). Accu- 
rate SEDs are computed for nine sources combining existing multi- 
wavelength surveys with new sub- mm data (Rengel et al. 2004b), 
and fitting a modified blackbody curve to the fluxes (Froebrich et 
al. 2003). We derived the bolometric temperature and luminosity, 
size of the envelope and sub-mm slope for each source. Furthermore, 
dust temperature distribution, SED, and an intensity map are de- 
rived self-consistently by the MC3D code. Physical parameters of 
the sources (e.g. envelope masses, density distributions, sizes and 
sublimation radius) were derived by finding the consistency between 
observed and modelled SEDs. 

We created a bolometric luminosity-temperature diagram for 
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Figure 2: Simulated and observed radial profiles as a function of the 
impact parameter b(AU) (left, short dashed and solid lines, respec- 
tively). The dotted line represents the la noise level at the edge of 
the map. Power-law index as function of time (right). 

the nine Class sources and compare the positions of the objects 
with an evolutionary model in order to derive the age and mass of 
the envelope of each source. 



4. Results 

In Fig. 1, the calculated temperature profile (left), and the ob- 
served and best simulated SEDs for L1448 C are shown as an exam- 
ple. Table 1 lists the characteristics of the best-fitting model pa- 
rameters. The envelope emission is calculated as an image from the 
MC3D code, and convolved by the SCUBA beam. The simulated 
intensity profile is calculated at 850 /mi and compared with the ob- 
served profile to test the symmetric model. 

Is there a correlation of the power-law index p with time? We 
plot the estimated values of p as a function of t. Fig. 2 suggests a 
density structure of p oc r~ 2 at younger ages, evolving to p oc r~ 3//2 
at later times. In order to investigate this with a larger spread, a 
Class 1 sample is mandatory. 



5. Conclusions 

Detailed physical interpretation of the internal structure of nine 
objects was carried out with intensive computer modeling. A simple 
spherically symmetric model envelope, and assumptions about den- 
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sity and dust distributions following the standard envelope model 
reproduce reasonably well the observed SED and the radial profiles 
of the sources. T oc r~ - 4 is a good approximation for the sample. 
The radial temperature distribution as function of distance, how- 
ever, departs significantly from the optically-thin assumption, an 
observational derivation of a single-power law q of 0.4 for radii < 10 
AU. Modelling results indicate a density profile well described by a 
power-law between p=1.5 to 2, which is expected by all of the co- 
llapse models and numerical studies. A density structure of p oc r~ 2 
at younger ages, evolving to p oc r~ 3 / 2 at later times was found. 

SED and radial profile fits can constrain physical parameters of 
the sources such as envelope masses, density distributions, sizes and 
sublimation radii. Nevertheless, observations at 10-300 pm, and the 
inclusion of outflow and disk, and other geometries will decrease the 
differences. 

Table 1: Best envelope fit results for the embedded sources. The 
temperature of the central object T* is set to 3500 K. x 2 quantifies 
the agreement between model and data. The age is estimated 
according to the model of Smith (2000) and given in 10 3 yrs. 



Object 


Tbol [K] 


L* [L ] 


Rsub [AU] 


P 


Rout [AU] 


age 


M e „„ [M ] 


x 2 


L 1448 NW 


27 


6 


5 


2.0 


4000 


13.8 


2.8 


6.4 


L 1448 C 


40 


11 


5 


1.6 


3000 


19.5 


1.7 


2.6 


RNO 15 FIR 


45 


8 


3 


1.6 


3500 


20.7 


1.0 


1.8 


NGC 1333 iras 1 


52 


13 


3 


1.5 


4500 


21.9 


1.5 


8.6 


NGC 1333 iras 2 


48 


68 


3 


1.5 


7000 


26.1 


3.0 


0.1 


HH211-mm 


30 


5 


3 


1.9 


6000 


13.8 


4.1 


4.5 


L1634 


42 


19 


3 


1.6 


6000 


21.9 


2.8 


1.1 


L 1641 N 


45 


85 


3 


1.5 


8500 


25.2 


7.0 


1.0 


L 1641 sms III 


50 


80 


3 


1.5 


10000 


27.6 


6.3 


0.6 
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